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Abstract 

High frequency stochastic resonance (SR) phenomena, associated with fluctuational 
transitions between coexisting periodic attractors, have been investigated experimentally 
in an electronic model of a single-well Duffing oscillator bistable in a nearly resonant field 
of frequency Up. It is shown that, with increasing noise intensity, the signal/noise ratio 
(SNR) for a signal due to a weak trial force of frequency Q ~ up at first decreases, then 
increases, and finally decreases again at higher noise intensities: behaviour similar to that 
observed previously for conventional (low frequency) SR in systems with static bistable 
potentials. The stochastic enhancement of the SNR of an additional signal at the mirror- 
reflected frequency 1^ — 2^1 is also observed, in accordance with theoretical predictions. 
Relationships with phenomena in nonlinear optics are discussed. 
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The first decade of research on stochastic resonance [1] (SR), in which the signal due to 
a weak (trial) periodic force in a nonlinear system can be optimally amplified by the intro- 
duction of external noise of appropriate intensity, has concentrated almost exclusively on 
systems with coexisting attractors corresponding to the minima of a symmetrical, static, 
bistable potential. In such cases, the noise-induced amplification arises through the oc- 
currence of fluctuational transitions between the attractors. For suitably chosen noise 
intensity, these become nearly periodic at the frequency of the trial force, with an average 
amplitude that can approach the half-separation of the attractors. This is the mechanism 
responsible for SR phenomena studied in connection with ice-ages, ring-lasers, electronic 
circuits, passive optical systems, electron spin resonance, sensory neurons, a magnetoe- 
lastic ribbon, and a laser with saturable absorber; we shall refer to it as conventional 
SR. 

More recently, efforts have been initiated to seek SR in other classes of systems having 
quite different kinds of attractors: it has now been identified in underdamped nonlinear 
oscillators that have single static attractors [2], in a bistable, chaotic, electrical circuit with 
two coexisting attractors of which one is a limit cycle and the other is chaotic [3] , and in a 
system with two coexisting stable limit cycles that have the same period and correspond to 
periodically forced vibrations of a damped nonlinear oscillator [4] ; the onset of a resonant 
absorption that increased extremely sharply with noise intensity - a phenomenon that 
has much in common with SR - had earlier been demonstrated theoretically [5] for the 
latter system. In this Letter we shall show that SR of this latter kind, while similar to 
conventional SR in some respects, also has a number of interesting features that distinguish 
it from earlier manifestations of the phenomenon and which have implications for four- 
wave mixing in nonlinear optics. We expect the ideas to be applicable to a large class 
of passive optically bistable systems and, in particular, to optically bistable microcavities 
[6]. 

The system that we consider is the nearly-resonantly-driven, underdamped, single- well 
Duffing oscillator with additive noise, 

q + 2Tq + cu*q + 7 g 3 = F cos(u F t) + /(f) (1) 
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r, \5u\ < uj f , 



^5uu > 0, 



5co = uof — 



(f(t)) = 0, (f(t)f(t>))=4TT5(t-t>) 

This system is of topical interest on account of its importance in nonlinear optics [7] 
and its relevance to experiments on a confined relativistic electron excited by cyclotron 
resonant radiation [8]. Provided that F 2 <C uOq(5uo 2 + r 2 )/|7|, and that the noise is 
weak, the resultant comparatively small amplitude [<C (^o/M) 5 ] oscillations of q(t) can 
conveniently be discussed in terms of the dimensionless parameters [5] 

" = r/|fo| - *=3a|&' ° = 8l?ir/&4r (2) 

which characterise, respectively, the frequency detuning, the strength of the main periodic 
field, and the noise intensity. The bistability [9] in which we are interested arises for a 
restricted range of rj and j3: within the triangular region bounded by the full lines of Fig 
1. Thus, as the amplitude of the periodic force is gradually increased from a small value 
at fixed frequency (see vertical line a - a'), the system moves from monostability (one 
small limit cycle), to bistability (two possible limit cycles of different radii), and then 
back again to monostability (one large limit cycle). The effect [5] of weak noise f(t) is 
to cause small vibrations about the attractors, and to induce occasional transitions (cf 
[10]) between them when the system is within the bistable regime. We shall see that 
SR phenomena occur in the close vicinity of the kinetic phase transition (KPT) line [11], 
shown dashed in Fig 1, where the populations of the two attractors are equal. 

Our principal aim is to consider the response of the system (1) to an additional weak 
trial force Acos(ilt + 0). The combined effects of dissipation and noise result in a steady 
statistical distribution, and the response of the system can therefore be described, in 
terms of linear response theory, by a susceptibility. The trial force beats with the main 
periodic force and thus gives rise to vibrations of the system, not only at Q, but also 
at the combination frequencies \Q ± 2uuf\ (and also at \Q ± (2n + 1)u>f\ in the case of 
nonlinearity of a general type). We are interested in the case where the strong and trial 
forces are both nearly resonant: that is, ujp and Q both lie close to uj . This is the case for 
which the response to the trial force is strongest. It is at its most pronounced at frequency 
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Q and at the nearest resonant combination, which for (1) is \Q — 2u>f\- The amplitudes 
of vibrations at these frequencies can be described respectively by susceptibilities x(^)> 
X(Q), so that the trial-force-induced modification of the coordinate q, averaged over noise, 
can be sought in the form 

5(q(t)) = ARe{ X (tt) exp[-ittt - i<f>] + X(Q) exp[i(2w F - Q)t] - i<f>]} (3) 

Within the KPT range, for Q close to cup, |Inr\;(f2)| becomes large and strongly noise- 
dependent [5]. The rapid rise in susceptibility with noise intensity corresponds precisely 
to SR since, according to (3), the squared amplitudes of the signals at frequencies fl and 
|Q — 2u>f\ (and the integrated powers of the corresponding peaks in the power spectrum) 
are 

S(Q) = \A*\xm 2 , S(\Q-2u F \) = \A"\X(Q)\" (4) 

An intuitive understanding of the mechanism of stochastic amplification can be gained 
by noting that the trial force modulates the driving force (and the coordinate q(t)) at 
frequency \Q — ujp\ and that, when \Q — cuf\ is small, the system responds almost adia- 
batically. In terms of the phase diagram Fig 1, the beat envelope then results in a slow 
vertical oscillation of the operating point. If the operating point is taken to be p, and its 
range of modulation p' — p" is set to straddle the KPT line as shown, and the noise in- 
tensity is optimally chosen, then the system will have a tendency to make inter-attractor 
transitions coherently, once per half-cycle of the beat frequency. The net effect of the 
noise is therefore to increase the modulation depth of the beat envelope of the response, 
thereby increasing the components of the signal at frequencies fl and \Q — 2u>p\- 

We are investigating the response of the system (1), and the variation of the sig- 
nal/noise ratio with a, through analogue experiments on an electronic model of (1), of 
which the relevant technical details were given in a recent conference report [4]. In terms 
of scaled units the circuit parameters were set, typically, to: 2r = 0.0397; lj = 1.00; 7 = 
0.1; uj f = 1.07200; Q = 1.07097; and, to seek SR near the KPT, F= 0.068 and the ampli- 
tude of the trial force A = 0.006. A spectral density of fluctuations of the coordinate q(t) 
(about (q(t)) for A=0)recorded with the above parameter values for a = 0.061 and 16384 
samples in each realization, is shown in Fig 2. The smooth background is the supernar- 
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row peak of [11], here broadened by noise (although its width remains very much smaller 
than 2r); delta function spikes, indicated by raised points [12] of the discrete spectrum, 
are clearly visible, not only at the trial force frequency (Q), but also as predicted at the 
mirror-reflected frequency (2up — Q). 

The signal/noise ratios R, determined in the usual way [1] from measurements of the 
delta spikes and the smooth background, are plotted (data points) as functions of noise 
intensity a in Fig 3 for f3 = 0.814, 77 = 0.236. It is immediately evident that there is 
a range of a within which R increases with a. It is also apparent that, for both the 
main and the mirror-reflected signals, the form of R(a) in Fig 3 is remarkably similar 
to that observed earlier [1] in the case of conventional SR. That is, R initially decreases 
with a, on account of the increase in its denominator. With further increase of a, the 
inter-attractor transitions come into play and become phase-coherent with the trial force 
to a high probability, giving rise to an increase in R through the stochastic amplification 
mechanism discussed above. Finally, for still larger a, R decreases again partly owing to 
the continuing rise in its denominator and partly because transitions are then occurring 
very frequently, within individual periods of the trial force, with a corresponding reduction 
in the proportion of the phase-coherent jumps that are responsible for the amplification. 

A quantitative theory of the phenomenon is readily constructed through an extension 
[13] of [5]. It leads to contributions to the susceptibilities from inter-attractor transitions 
of the form 

(5) 

where Wi,w 2 are the populations of the attractors 1, 2 and W r 12 , W 2 \ are the probabilities 
of transitions between them, which are of the activation type Wij oc exp(— Ri/a). (The 
Ui,u*, which define the positions of the attractors in the rotating coordinate frame, can 
be regarded as constants for given 77, ($). It is evident that the contributions (5) come 
into play if, and only if, the system is within the KPT range where the populations of the 
unperturbed attractors are comparable: otherwise, the factor w±w 2 oc exp(— \R 1 — R 2 \/a) 
will be exponentially small. Within the KPT range, however, the susceptibilities will be 
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large because they are proportional to the large factor \fii — /U 2 |/o;; the rapid increase 
of Wij with noise intensity then ensures that there will be a range of a in which both 
susceptibilities increase very rapidly with a, consistent with the experiments. The full 
theory [13], including the effect of intra-attractor vibrations, leads to the curves of Fig 
3. Given the large systematic errors inherent in the measurements - arising e.g. from Su 
(1), a small difference between large quantities which, in (3 (2), is then raised to its third 
power - the agreement between theory and experiment can be considered excellent. 

In conclusion, we would emphasize, first, that, in contrast to earlier forms of bistable 
SR [1, 3], stochastic amplification occurs here not at the relatively low frequency of the 
quasi-periodic inter-attractor hopping but, rather, at Q close to the much higher (tunable) 
frequency up of the main periodic driving force. To emphasize the distinction, it seems 
appropriate to refer to the new form of SR as high frequency stochastic resonance (HFSR). 

Secondly, we draw attention to the relationship of HFSR to four-wave mixing in non- 
linear optics [14]. In effect, our results correspond to noise-enhanced amplification of 
the signal wave, and noise-enhanced generation of the idler wave. The mechanisms are 
resonant and, although they have the appearance of four-wave mixing, they actually cor- 
respond to multiple-wave processes: in terms of quantum optics, the oscillator absorbs 
and re-emits many quanta of the strong field. The very high amplification/generation 
coefficients arise partly from their resonant character and partly from the fact that the 
signal sizes correspond, not to the amplitudes of vibrations about the attractors but, as 
usual in bistable SR, to the "distance" between the attractors (to their coordinate sep- 
aration for conventional SR, and approximately to their difference in amplitude in the 
present case). 

Finally, our prediction and demonstration of HFSR for periodic attractors, and its 
similarity (Fig 3) to conventional SR, leads to a broader and more general perception of the 
physical nature of bistable SR. Like the onset of supernarrow peaks in the power spectra, 
conventional SR [1] and HFSR are both critical phenomena that arise in the range of the 
KPT. While HFSR is to be anticipated for coexisting stable limit cycles with equal periods, 
low-frequency SR is a more robust effect. It arises, not only for systems fluctuating in 
simple double-well potentials but also for systems where one (or both) of the attractors is 
chaotic [3]; we may infer from [5] that, in general, low frequency SR is also to be anticipated 
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for periodic attractors [although not for (1), where the centres of the forced vibrations 
are independent of amplitude]. Since noise gives rise to fluctuational hopping between 
any type(s) of attractors, it seems reasonable to conclude that SR is actually a quite 
general phenomenon characteristic of all systems with coexisting attractors, regardless of 
the nature of those attractors, provided only that the system lies within its KPT range. 

This research has been supported by the Science and Engineering Research Council 
(United Kingdom), by the Royal Society of London and by the European Community. 
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Figure Captions 



1. Phase diagram for (1) in terms of reduced parameters (2): the cuts a-a', p'-p-p" are 
discussed in the text. 

2. Power spectral spectral density Q(uS) measured for the electronic model, with the 
contents of each FFT memory address shown as a separate data point on a highly 
expanded abscissa; a smooth curve, peaking at up, is drawn through the background 
spectrum; vertical lines indicate the delta spikes resulting from the trial force. 

3. The signal/noise ratio R of the response of the system (1) to a weak trial force at 
frequency Q, as a function of noise intensity a, in experiment and theory: at the 
trial frequency Q (circle data and associated theoretical curve); and at the "mirror- 
reflected" frequency (2uf — ^) (squares). For noise intensities near those of the 
maxima in R(a), the asymptotic theory is only qualitative and so the curves are 
shown dotted. 
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